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Abstract. We present decompositions of the rotation curves of three spiral galaxies at redshifts z ~ 0.7 and 1 into contributions 
by their bulges, disks, and dark halos, respectively. In order to set constraints on the degeneracy of the decompositions we inter- 
pret the morphology of the spiral structures quantitatively in the framework of density wave theory. Galaxy models constrained 
in such a way show that the distant galaxies, which are much younger than nearby galaxies, have very likely 'maximum disks', 
i.e. are dominated in their inner parts by baryonic matter. We argue that current theories of the cosmogony of galaxies must 
allow for these types of galaxies. 
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1. Introduction 

The rotation curves of spiral galaxies, in particular the HI - 
rotation curves, have incontrovertibly revealed that the outer 
parts of the galaxies are dominated by dark matter. For in- 
stance, van Albada et al. (1985) have measured the rotation 
curve of NGC3198 in HI extending radially outwards to 11 
disk scale lengths, much beyond the optical disk. At such a 
galactocentric radius the rotation curve of a single disk would 
have fallen to half of its peak value at 2.2 radial scale lengths 
(cf. equation 5), whereas the observed rotation curve stays flat. 
However, there is an intense debate in the literature if and to 
what degree the inner parts of the galaxies, where their op- 
tical disks reside, are dominated by dark matter as well. The 
various arguments and the present state of this so called 'max- 
imum disk' problem have been reviewed recently by Bosma 
(2004). The reason for this dispute is the notorious degeneracy 
of the decomposition of the rotation curves into the contribu- 
tions by the various constituents of the galaxies such as bulge, 
stellar disk, interstellar gas disk, and dark halo, respectively. 
Thus, further constraints on the decomposition of the rotation 
curves are needed. Using arguments of the density wave theory 
of galactic spiral arms Fuchs (2003) has pointed out that the im- 
plied dynamics of the disks might provide such constraints for 
modelling spiral galaxies. Its application to the nearby bright 



Send offprint requests to: B. Fuchs 
* e-mail:fuchs @ ari.uni-heidelberg.de 



galaxies NGC2985 and NGC3198 points towards maximum 
disk models in both cases, i.e. the contributions of the disks to 
the rotation curves are at the maximum allowed by the data. 
This would imply that baryons dominate completely the inner 
parts of the galaxies. 

If this constraint is applied to the decomposition of the 
rotation curves of redshifted galaxies, this allows insight into 
the cosmogony of spiral galaxies. For this purpose Fuchs, 
Mollenhoff & Heidt (1998) have analyzed the optical rotation 
curves of two spiral galaxies at redshifts z — 0.15 and z = 0.48, 
respectively, which had been imaged with HST and showed 
clearly discernible spiral structures. The rotation curves were 
measured with the Keck telescope (Vogt et al. 1996). Fuchs et 
al. (1998) concluded that the disks of these galaxies are sub- 
maximal. This has to be modified, though, because they as- 
sumed flat rotation curves for the galaxies. Since the rotation 
curves span radially only two and three radial scale lengths, re- 
spectively, the rotation curves are still rising over these ranges, 
which had not been taken into account. The disks were mod- 
elled as exponential disks, and their rotation curves rise until 
2.2 radial scale lengths and develop only then a flat plateau 
(cf. equation 5). This affects the density wave theory argument 
through the coefficient X(A) in equation (8) below, and, when 
properly taken into account, the disks of the two redshifted 
galaxies seem to be maximum disks as well (Leeuwin & Fuchs, 
priv. com.). 
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Fig. 1. Upper panel: HST images of FDF2484 (z = 0.65), FDF2174 (z = 0.68), and FDF4922 (z = 0.97). North is up and east 
is left. The scale is 0.05 arcsec/pix, the image size is 6" x 6". 

Lower panel: Corresponding residual images after subtraction of 2-dim models for disk and bulge. The remaining structures 
show the dominant bisymmetric arms. White, grey, and black colors corresponds to positive, zero, and negative residuals, respec- 
tively 



In the present paper we investigate spiral galaxies identi- 
fied in the FORS Deep Field (hereafter referred to as FDF, 
Appenzeller et al. 2000). The photometric survey is described 
in detail by Heidt et al. (2003). Ziegler et al. (2002) and Bohm 
et al. (2004) observed the kinematics of a large sample of FDF- 
galaxies in the redshift range z = 0.1 - 1 . Three of the galaxies 
at z > 0.5 turn out to be suitable for a quantitative analysis of 
their morphologies and kinematics and allow to extend our pre- 
vious studies out to redshifts z ~ 1- The data are described in 
section 2. In sections 3 and 4 we construct dynamical models 
for the galaxies. As discussed in the final section the dynam- 
ical modelling seems to indicate maximum disks also in field 
galaxies in that redshift range. 

2. Observations 

Bohm et al. (2004) have measured spatially resolved rotation 
curves of 77 galaxies in the FDF using the FORS spectro- 
graph at VLT. 36 of these rotation curves were classified as high 
quality in terms of large radial extent and small asymmetries. 
Moreover, Ziegler and collaborators (Bohm et al., in prepara- 
tion) have obtained HST images of the FORS Deep Field with 



Table 1. Redshifts, absolute rest frame magnitudes, 

seeing (FWHM), maximum rotation velocities, 

and scale conversion factors of the observed galaxies 



FDF 


redshift 




seeing 




scale 






mag 


arcsec 


km/s 


kpc/arcsec 


2484 


0.6535 


-21.07 


0.80 


169 


6.14 


2174 


0.6798 


-21.00 


0.51 


175 


6.26 


4922 


0.9731 


-21.78 


0.43 


156 


7.20 



the ACS wide field camera. The ~ 6 x 6 arcmin field-of-view 
of the FDF was observed within 4 visits. Each visit was splitted 
into two exposures of 1 1 80 s integration time using the F8 14W 
filter. We have examined the ACS image of each high quality 
object in the redshift range z = 0.5 - 1 and found three galax- 
ies showing discernible spiral structure: FDF 2484 at redshift 
z ~ 0.65, FDF 2174 at z ~ 0.68, and FDF 4922 at z ~ 0.97, 
respectively. The HST images of the galaxies are reproduced 
in the upper panels of Fig. [2 As can be seen from the residual 
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Table 2. Photometric parameters from the bulge-disk decom- 
positions. Zbo and Sjo are the projected central flux-densities 
of bulge and disk, measured in analog-digital units (adu) per 
pixel. /?b and R d are the characteristic scale lengths, i is the in- 
clination of the disk and 8 is the angle between spectrograph 
slit and apparent major axis of the galaxy. 



FDF 










i 


e 




adu/pix 


adu/pix 


arc sec 


arc sec 


deg 


deg 


2484 


157 


133 


0.13 


0.65 


54 


-12 


2174 


182 


251 


0.12 


0.50 


67 


-4 


4922 


64 


62 


0.22 


0.82 


71 


14 



images in the lower panels all galaxies show well developed 
two-armed spirals. The m — 2 arm mode was confirmed by 
a Fourier analysis of the arm pattern as described in the next 
section (Fig.|2ji. 

The optical rotation curves are shown as data points in 
Figs.|5]to|5]in section 4. The seeing during the spectroscopic 
observations is summarized in Table 1 . For the determination 
of the misalignment between the galaxies' apparent major axes 
and the slit directions, we took into account a slight rotation 
of the HST/ACS frames with respect to the FORS spectra. This 
tilt has been derived by coordinate transformations between the 
two reference systems and amounts to approximately -4 deg. 
The last column in Table 1 gives scale conversion factors from 
arcsec to kpc. These have been determined in the same way 
as by Vogt et al. (1996) adopting a cosmological model with 
qo = 0, which at the redshifts discussed here is similar to the 
concordance model with £2 m = 0.3 and D.\ = 0.7. We assumed 
a Hubble constant of Ho = 75 km s _1 Mpc -1 . 

As can be seen from Figs. 3 to 5 the rotation curves have 
been measured with good accuracy. The methods to extract the 
rotation curves from the spectra and to estimate the errors are 
described in detail in Bohm et al. (2004). The rotation curves 
are very regular and quite symmetric. Taken together with the 
unperturbed optical appearances of the galaxies this should 
allow in our view a reliable comparison with model rotation 
curves in order to draw inferences about the dark matter con- 
tent. 



3. Bulge-disk decomposition 

The surface-brightness distributions of the galaxies were mod- 
elled in the same way as in Fuchs, Mollenhoff & Heidt (1998). 
A thin disk was assumed which has elliptical isophotes when 
seen under an inclination angle i. The radial brightness profiles 
were fitted by an exponential function 



Sd(/?) = Sd exp(-/?//? d ) 



(1) 



with a scale length R&. For reasons of the easier handling of the 
dynamical models we adopt instead of a de Vaucouleurs law 



for the surface brightness profiles of the bulges softened power 
laws 



SbW = Em 



1 + ~7 



-2.5/2 



(2) 



where R^, denotes the core radius. The corresponding three- 
dimensional distribution has the same form as equation (2), 
but with an exponent lowered to -3.5/2. The central density 
of the three-dimensional distribution pbo is given by pbo = 
0.572 Zao/^b- The models (1) and (2) were fitted simulta- 
neously to the data using the two-dimensional fitting code 
of Mollenhoff (2004). The parameters Eb0>£d0>^b>^d> h and 
the position angle of the major axis of the disk are free fit- 
parameters. The photometric parameters of the best fits are 
given in Table 2. The angle 9 in Table 2 is the position angle of 
the slit relative to the major axis of the inclined galaxy. 

The fit to the surface-brightness distribution describes the 
axisymmetric structure of disk and bulge. Therefore the spi- 
ral arms appear clearly in the residual image = galaxy - model 
(Fig-E lower panels). The residual images were transformed to 
face-on view and the spiral arm pattern was Fourier-analyzed. 
Fig. |2 shows the radial variation of the coefficients m-2, 3, 4 
in the galaxies. The two-armed pattern m-2 is dominant in all 
cases. The image of FDF2174 shows that the spiral extends ra- 
dially out to about 1.3 arcsecs. This corresponds in the diagram 
of the m-2 Fourier coefficients in Fig.|2]to the indentation at 
a radius of about 25 pixel lengths. The very large amplitude 
of the m-2 Fourier coefficients beyond that radius is caused 
by two outer patches lying nearly diagonally towards northeast 
and southwest, which are, however in our view, not related to 
the spiral structure (cf. Fig.Q. But as evidenced by the rotation 
curve, these outer regions do belong to the optical disk of the 
galaxy. 

4. Decomposition of the rotation curves 

We adopt model rotation curves of the form 

= v lb( R ) + <dW + »&(*) . (3) 

where v c w c ,d, an d v c ^ denote the contributions by the bulge, 
disk, and dark halo, respectively. The bulge contribution is 
given by 



f 

Jo 



drr 1 



r 2 \-3-5/2 

1 + - 



(4) 



where G denotes the constant of gravitation. The rotation 
curve of an infinitesimally thin exponential disk is given by 
Freeman's (1970) formula, 



v 2 c /R) = 4nGl d0 R d x 2 ( I (x)K (x) - h(x)K x (x) ) , 



(5) 



where x is an abbreviation for x = R/2R d and / and K de- 
note Bessel functions. In equations (4) and (5) we use the scale 
lengths determined in the previous section and the ratios of pbo 
fLdo found by the quantitative modelling of the surface bright- 
ness photometry. Only the mass-to-light ratio is a free parame- 
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Fig. 2. Fourier coefficients (in arbitrary units) m = 2 (solid line), m = 3 (dashed line), and m = 4 (dotted line) of the spiral arm 
patterns in FDF2484, FDF2174, and FDF 4922, respectively. The coefficient m = 2 is dominant in all cases. 



ter which is determined as explained below. The dark halos are 
modelled as quasi-isothermal spheres 



Ph(R) = PhO 



r 2 

1 + K 



(6) 



which leads to a contribution to the rotation curve of the form 



v 2 c d (R) = 47rG Ph0 R 2 h (l - y alCtan { R \)) ■ 



(7) 



Such density profiles of dark halos with homogeneous cores 
fit usually the rotation curves of galaxies better (de Blok, 
McGaugh & Rubin 2001, de Blok & Bosma 2002, Gentile et 
al. 2004) than centrally cusped profiles as proposed, for in- 
stance, by Navarro, Frenk & White (1997). The velocity fields 
described by the rotation curve model (3) are projected onto 
the sky adopting the inclination angles determined in the previ- 
ous section. Next, the line-of-sight components of the rotation 
velocities are calculated and weighted by the surface bright- 
ness. The radial velocity fields were then convolved with two- 
dimensional Gaussians in order to model the blurring of the 
velocity fields by seeing. Finally the velocity field models were 
masked according to a slit width of 1 arcsec and the position 
angles of the slit given in Table 2 (cf. also Fig. 4 of Bohm et al. 
2004). The synthetic rotation curves determined this way can 
be directly compared with the observed rotation curves. 

Since the decomposition of the rotation curves is com- 
pletely degenerate, we set constraints on the decomposition 
using arguments of the density wave theory of galactic spi- 
ral arms. The density wave theory predicts that spiral den- 
sity waves grow preferentially with azimuthal wavelengths 
(Toomre 1981, Fuchs 2001) 



A = X(A)A clit = X(A) 



47!- 2 G£ d 



(8) 



where k denotes the epicyclic frequency of the stellar orbits, 
k = y/2(v c /R) VI + din v c /dlnR. The coefficient X depends 
on the slope of the rotation curve measured by Oort's constant 
A (Fuchs 2001). For a flat rotation curve its value is X = 2. The 



Table 3. Dynamical parameters from the decomposition of the 
rotation curves, pbo is the central spatial density of the bulge 
(in Alo/pc 3 ). E d is the central surface density of the disk 
(M /pc 2 ), and pho is the central spatial density of the halo. The 
next 3 columns give the mass of bulge, disk, and halo, within 
radii 12, 10, and 15 kpc for FDF 2484, 2174, and 4922, resp. 
The last column is the resulting rest frame B mass-to-light ratio 
calculated using total disk mass estimates. 



FDF 


PbO 

Mo 
pel 


2d 

Mo 
pc 2 


PhO 

Mp 
pf? 


Mb 


Mi 
\0 9 M o 


Mb 


Al 
X 

Me 
Xbo 


2484 


0.41 


480 




4.6 


39 




1.5 




0.38 


450 


0.0016 


4.6 


36 


6.5 


1.4 


2174 


0.38 


700 




3.7 


36 




1.4 




0.33 


600 


0.002 


3.7 


31 


5.4 


1.2 


4922 


0.13 


350 




11 


55 




1.3 




0.12 


310 


0.0013 


11 


49 


8.4 


1.2 



expected number of spiral arms m is obviously given by how 
often the circumferential wavelength A fits onto the annulus, 



2nR 



(9) 



As was pointed out first by Athanassoula (1988) and 
Athanassoula et al. (1987) equations (9) and (8) can be used to 
determine the surface density of galactic disks once the rotation 
curve v c (R) and thus the epicyclic frequency k are known. In 
this way the morphological appearance of spiral galaxies con- 
strains the decomposition of the rotation curve into its various 
contributions. 

We have tested the prediction of equation (9) on the numer- 
ical simulation of the dynamical evolution of a galactic disk by 
Sparke & Sellwood (1987). The mass of the disk and the ro- 
tation curve are known from the setup of the simulation and 
we find that equation (9) predicts indeed a two-armed spiral as 
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Fig. 3. Upper panel: Model rotation curves of FDF 2484 fit- 
ted to the radial velocity data of Bohm et al. (2004). The solid 
line is the model without a dark halo, the dotted line shows the 
corresponding rotation curve when a dark halo component was 
included. The rotation curves are practically identical. 
Lower panel: Expected number of spiral arms in the radial 
range where spiral structure can be seen. The dotted curve is 
again for the model with a dark halo included. 
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Fig. 4. Same as Fig.|3] but for FDF2174. 

seen in the first snapshots of the simulation before it develops 
a bar. 

5. Results and discussion 

In the upper panels of Figs. [3]to[5]we present model rotation 
curves fitted to the observed rotation curves of the three galax- 
ies assuming no dark halos. As can be seen from the figures 
the models fit well to the data. The corresponding dynamical 
parameters are listed in Table 3. The expected number of spi- 
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Fig. 5. Same as Fig.|3 but for FDF 4922. 

ral arms estimated using equation (9) are shown as solid lines 
in the lower panels of Figs.[3]to[5] In each case the prediction 
is a two-armed spiral, exactly as observed. This hardly allows 
massive inner dark halos. Nevertheless, we also performed fits 
with dark halos. 

Unfortunately, the rotation curves span radially only be- 
tween two and four disk scale-lengths. The core radii of the 
dark halos are not well constrained in such cases (van der Kruit 
1995). We have therefore adopted rather arbitrarily for each 
galaxy a core radius of 10 kpc. The central halo densities were 
chosen so that the expected number of spiral arms is still in 
accordance with the observed morphologies of the galaxies, 
i.e. two-armed spirals. The resulting fits are shown as dotted 
lines in Figs.[3]to|5] As can be seen from the figures they are of 
the same quality as in the case of the bulge-disk models with- 
out dark matter. The dark halo contributions could be raised 
and simultaneously the disk contributions lowered further, and 
would lead to equally good fits to the rotation curves. However, 
the expected number of spiral arms would then increase ac- 
cording to equation (9) as m oc S d ' and would fit no longer to 
the observed morphologies of the galaxies. Thus, applying the 
density wave theory constraint to the decomposition of the ro- 
tation curves the disks of the three galaxies in the z = 0.7 to 1 
redshift range we find that the galaxies seem to have maximum 
disks similar to the nearby galaxies. This does not mean that we 
argue against the presence of dark matter in the galaxies. The 
observed rotation curve of FDF 2 174, for example, spans al- 
ready about four radial scale lengths and shows no indications 
of falling at the outer radii. This would be the case if there were 
only a single disk. Instead the flat shape of the rotation curve 
implies according to our model a ratio of A / lh/(A / lb + Aid) 
0.24 at galactocentric radius R = 12.5 kpc. 

Using the absolute rest frame luminosities of the FDF 
galaxies derived by Bohm et al. (2004), transformed to the 
value of Ho adopted here, we have determined the mass-to- 
light ratios shown in Table 3. The dynamically determined 
mass-to-light ratios are in good agreement with stellar mass- 
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to-light ratio estimates based on population synthesis models 
for galaxies in this this redshift range (Dickinson et al. 2003, 
Drory et al. 2004). 

We have demonstrated that there exist galaxies out at red- 
shifts of z » 1 which are dominated by baryons in their inner 
parts. Obviously this must reflect on theories of the cosmogony 
of the galaxies. It would be interesting if the present paradigm 
of A cold dark mater cosmology can account for such galaxies. 
For instance, the recent state-of-the-art simulation of galaxy 
formation by Abadi et al. (2003) lead to a galaxy model domi- 
nated by dark matter through all phases of its evolution, which 
is at variance to our findings. Admittedly the selection criterion 
of visible spiral structure might introduce a bias in the galaxy 
sample which we have analyzed. Since the galaxies have orga- 
nized themselves quite early by self gravity, they have proba- 
bly not experienced a major merger event and could settle into 
equilibrium. 
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